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O- ! Abstract 

o 

■ We apply a microcanonical statistical model to investigate hadron produc- 
tion in pp collisions. The parameters of the model are the energy E and the 



i-O \ volume V of the system, which we determine via fitting the average multi- 

plicity of charged pions, protons and antiprotons in pp collisions at different 
Oh. collision energies. We then make predictions of mean multiplicities and mean 

transverse momenta of all identified hadrons. Our predictions on nonstrange 
hadrons are in good agreement with the data, the mean transverse momenta of 
strange hadron as well. However, the mean multiplicities of strange hadrons 
are overpredicted. This agrees with canonical and grandcanonical studies, 
where a strange suppression factor is needed. We also investigate the influ- 
ence of event-by-event fluctuations of the E parameter. 



1 Introduction 

In pp collisions at high energies a multitude of hadrons is produced. In contradis- 
tinction to the pp collisions at low energies even effective theories are not able to 
provide the matrix elements for these reactions and therefore a calculation of the 
cross section is beyond the present possibilities of particle physics. In addition, even 
at moderate energies, many different particles and resonances may be created and 
therefore the number of different final states becomes huge. 

* Fellow of Alexander von Humboldt Foundation 
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In this situation, statistical approaches may be of great help [U El- It was 
Hagedorn who noticed that the transverse mass distributions in high energy hadron- 
hadron collisions show a common slope for all observed particles[3j. This may be 
interpreted as a strong hint that it is not the individual matrix elements but phase 
space who governs the reaction. Therefore Hagedorn introduced statistical methods 
into the strong interaction physics in order to calculate the momentum spectra of 
the produced particles and the production of strange particles. 

Later, after statistical models have been successfully applied to relativistic heavy 
ion collisions [4j El O IZl IHl IHl CHI [TTJ , Becattini and Heinz |J2] came back to the 
statistical description of elementary pp and pp reactions and used a canonical model 
(in which the multiplicity of hadrons M is a function of volume and temperature 
M(V,T)) in order to figure out whether the particle multiplicities predicted by this 
approach are in agreement with the (in the meantime very detailed) experimental 
results. For a center of mass energy of around 20 GeV they found for non strange 
particles a very good agreement between statistical model predictions and data 
assuming that the particles are produced by a hadronic fireball with a temperature 
of T = 170 MeV. The strange particles, however, escaped from this systematics being 
suppressed by factors of the order of two to five. Becattini and Heinz coped with 
this situation by introducing a 75 factor into the partition sum which was adjusted 
to reproduce best the multiplicity of strange particles as well. 

Statistical models are classified according to the implementation of conservation 
laws: 

• microcanonical: both, material conservation laws (Q, B, S, C, ■ ■ •) and mo- 
tional conservation laws (E, p , J , ■ ■ •), hold exactly. 

• canonical: material conservation laws hold exactly, but motional conservation 
laws hold on average (a temperature is introduced). 

• grand-canonical: both material conservation laws and motional conservation 
laws hold on the average (temperature and chemical potentials introduced). 

The intensive physical quantities such as particle density and average transverse 
momentum are independent of volume in the grand-canonical calculation , while 
they depend on volume in both canonical and microcanonical calculations. What 
one naively expects is that the microcanonical ensemble must be used for very small 
volumes, for intermediate volumes the canonical ensemble should be a good approxi- 
mation, while for very large volumes the grand-canonical ensemble can be employed. 
A numerical study of volume effects in paper[Hj tells us how big the volumes need 
to be in order to make the grand-canonical ensembles applicable. The comparison 
between the microcanonical and the canonical treatment in paperfHj shows a very 
good agreement in particle yields, when the same volume and energy density are 
used, and the strangeness suppression is canceled in the canonical calculation. 
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In this paper, first we ignore the fluctuations of microcanonical parameters and 
try to fix the microcanonical parameters, energy E and volume V, from fitting Air 
yields of protons, antiprotons and charged pions from pp collisions. The one-to-one 
relation between the collision energy \fs and a pair of microcanonical parameters 
E and V makes a link between the pp experiments and the microcanonical ap- 
proaches (or more generally, the statistical ensembles). One can easily judge if 
grand-canonical ensembles can describe pp collisions at any given energy; one can 
also transform the fitting results to the canonical case and find the corresponding 
temperature and volume of pp collisions at any energy. 

Then we study the effect from the fluctuations of the microcanonical energy 
parameter at a collision energy of 200 GeV, to check how reliable it is to fix micro- 
canonical parameters without energy fluctuations. 

Finally, we would like to make a comparison between statistical models and 
string models in describing pp collisions. This microcanonical model and this fitting 
work will provide us a bridge to compare the two classes of models and help us to 
understand the reaction dynamics. In principle, one can consider a string as an 
ensemble of fireballs, which may be considered as one effective fireball, when only 
total multiplicities are considered. 



2 The approach 

We consider the final state of a proton-proton collision as a "cluster" , "droplet" 
or "fireball" characterized by its volume V (the sum of individual proper volumes), 
its energy E (the sum of all the cluster masses) and the net flavour content Q = 
(N u — N u ,N d — Nj,N s — Ng), decaying "statistically" according to phase space. 
More precisely, the probability of a cluster to hadronize into a configuration K = 
{hi,pi; . . . ; h n ,p n } of hadrons hi with four momenta pi is given by the micro- 
canonical partition function Q(K) of an ideal, relativistic gas of the n hadrons 

CHI, 

\rn n i n 

\ / i=l a&S a ' i=l 

with £j = \Jmf + pf being the energy, and Pi the 3-momentum of particle i. The 
term ensures flavour conservation; is the flavour vector of hadron i. The 

symbol S represents the set of hadron species considered: we take S to contain 
the pseudoscalar and vector mesons (ir, K,r),r)' , p, K*,ui,<fi) and the lowest spin-| 
and spin-| baryons (N, A, S, H, A, £*, 5*, fl) and the corresponding antibaryons. n a 
is the number of hadrons of species a, and ^ is the degeneracy of particle i. We 
generate randomly configurations K according to the probability distribution fl(K). 
For the details see ref. [TH]. The Monte Carlo technique allows to calculate mean 
values of observables as 

K K' 



3 




Figure 1: The 4n multiplicities of 7r + ,7r + , proton, antiproton produced in a pp 
collision as a function of y/s. The full and dashed lines show the result of the \ 2 fit. 
The points are data from [T5] . 



where means summation over all possible configurations and integration over the 
Pi variables. A(K) is some observable assigned to each configuration, as for example 
the number M h (K) of hadrons of species h present in K. Since A depends on E and 
V, we usually write A(E, V). Q is not mentioned, since we only study pp scattering 
here, therefore Q is always (4,2,0). 

Let us consider the hadron multiplicity Mh(E,V). This quantity is used to 
determine the energy E(y/s) and the volume V(y/s) which reproduces best the 
measured multiplicity of some selected hadrons in pp collisions at a given y/s . This 
is achieved by minimizing x 2 ' 

where M exp j(y/s) and Oj are the experimentally measured multiplicity and its error 
of the particle species j in pp collisions at an energy of y/s. 

We start out our investigation by taking as input the most copiously produced 
particles (j = p,p, n + , The data have been taken from [TH]. Whenever the data 
are not available, the extrapolation of multiplicities by Antinucci [TH] is used. Fig.l 
displays the results of our fit procedure in comparison with the experimental data. 
We observe that these 4 particle species can be quite well described by a common 
value of E(y/s) and V(y/s). 

Fig. 2 shows E(y/s) and V(y/s), and fig. 3 the energy density e(y/s) = E (y/s) /V (y/s), 
which we obtain as the result of our fit. Both energy and volume increase with y/s 
but rather different, as the energy density shows. We parameterize the energy and 
volume dependence on the collision energy y/s in eq. (JTJ). 

E/GeV = -3.8 + 3.761nv^ + 6.4/v^ 
y/fm 3 = -30.0376 + 14.93^^-0.013^ (1) 

where y/s is in unit GeV. Below y/s = 8 GeV the fit produces volumes below 2fm 3 
which cannot be interpreted physically. Above y/s = 8 GeV the volume increases 
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Figure 2: The dependence of micro-canonical parameters E (left) and volume V 
(right) on the collision energy \fs . The parameterization described as eq. (pQ) is 
plotted as dashed lines. 




Figure 3: The dependence of energy density e = E /V on the collision energy \/s. 
The dashed line corresponds to constant energy density 0.342 GeV/fm 3 which comes 
from a canonical calculation jl2*| IT1]. 



very fast as compared to the energy giving rise to a decrease in the energy density 
until - around \fs = 200 GeV - the expected saturation sets in and the energy density 
becomes constant. In view of the large volume observed for these large energies the 
density of the different particles does not change anymore [14f and therefore the 
particle ratios stay constant above this energy. 

The quality of the fit can be judged from fig. 4 where we have plotted the x 2 
values obtained for different values of E and V and for \fs = 200 GeV. We see that 
the energy variation is quite small whereas the volume varies more. Nevertheless 
the energy density is rather well defined. 

After having fitted the E(y/s) and V(\/s) in using p,p,ii + and n~ data we can 
now use these fitted values to predict the multiplicity of other hadrons. This study 
we start in Fig. 5, where we present the multiplicity of 7T° and p°. For these particles 
experimental data are available. We see that the absolute value as well as the trend 
of the experimental data is quite well reproduced. The result for those hadrons, for 
which no or only few data are available, is displayed in fig. 6. As one can see the 
overall agreement is remarkable. We would like to mention that we have as well 
made a \ 2 fit using as input the measured multiplicities of p,p and p°. The results 
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p+p at \Js=200GeV, % 2 plot 
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Figure 4: The value of y 2 for different values of E and V for a pp reaction at 
y/s = 200 GeV. 



for E(y/s) and V(y/s) differ only marginally. 



3 Strange particles 

With the parameters E(y/s) and V(y/s) which we have obtained from the fit of the 
p, p, 7r + and 7r~ multiplicities, we can as well calculate the multiplicity of strange 
particles or particles with hidden strangeness. The results of these fits are presented 
in Fig. 7. As we can see immediately the results for those particles are not at all in 
agreement with the data. A and multiplicities are off by a factor of 3-5 roughly, for 




Figure 5: Prediction of the 7r°and p multiplicity in pp collisions as a function of 
y/s. The result of the calculation is compared to the data 



6 




Figure 6: Predictions of the multiplicities of non strange hadrons in pp collisions 
as a function of yfs. We have plotted, if available, also the data points for y/s = 
27.5 GeV [IZ| 
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Figure 7: Prediction of the multiplicity of strange hadrons in pp collisions as a 
function of y/s. The result of the calculation is compared to the data [TK] . 

the A the situation may be similar but the spread of the experimental data does not 
allow for a conclusion yet. Only the kaons come closer to the experimental values. 
Although at lower energies a part of this deviation may come from the fact that 
in our Monte Carlo procedure weak decays are neglected and therefore K° and K° 
are the particle states which are treated, at higher energies this is not of concern 
anymore. At y/s = 53 GeV, we find - as in experiment - that Kl = K$ and hence 
K° = K°. Therefore, as in experiment, one finds that the strangeness contained 
in A, A, K + and K~ adds up to zero. The absolute numbers are, however, rather 
different: experimentally one finds .41 K + , 0.29 K~ and 0.12 A (THl 1X0] , whereas 
the fit yields 1.16 K+, 0.66 K~ and 0.58 A. 

One is tempted to try to fit the strange particles separately. The result at 
large y/s is that, in contradiction to experiment, more K° then K° are produced. 
Consequently, the strangeness in A, A, K + and K~ does not add up to zero and the 
fit is far away from the data. Thus we have to conclude that the strange particle 
multiplicities cannot be described in a phase space approach using the parameters 
one obtains from the fit of non strange particles, and that there is no understanding 
presently why the suppression factor is rather different for the different hadrons. 

As mentioned above, in the past there has been proposed to use an additional 
parameter, j s , in order to describe the strangeness suppression. This parameter 
has been interpreted as a hint that the volume in which strangeness neutrality has 
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to be guaranteed is small as compared to the volume of the system. However, a 
detailed comparison of the multiplicity of all strange particles with the data shows 
[T2] that one additional parameter alone is not sufficient to describe the measured 
multiplicities of the different strange particles in a phase space approach to pp 
collisions. 

4 Transverse Momenta 

Phase space calculations predict not only particle multiplicities but also the mo- 
menta of the produced particles. The average transverse momenta of the produced 
particle give a good check whether the energy density obtained in the fit can really 
be interpreted as the energy density of a hadron gas. Fig. 8 shows the average 
transverse momenta in comparison with the experimental data ^H]- We see that 
over the whole range of beam energies the average transverse momenta are in good 
agreement with the data. This confirms that the partition of the available energy 
into energy for particle production and kinetic energy is correctly reproduced in the 
phase space calculation. It is remarkable that the average transverse momenta of 
strange particles is correctly predicted. 

5 Energy Fluctuations 

Up to now we have assumed that for a given center of mass energy, the energy of 
the droplet E has an unique value, given in fig. 2. This is of course not a realistic 
assumption. Most probably the energy varies from event to event but little is known 
about the form of this fluctuation. The only quantity for which data are available 
is the multiplicity distribution of charged particles, which has been the subject of 
an extensive discussion in the seventies due to the finding of a scaling law, called 
KNO scaling. Of course, one can try now to find an energy distribution which yields 
the experimental charge particle distribution but this relation is not unique and 
therefore little may be learnt. 

It has also been suggested to replace the microcanonical ensemble calculation, 
presented here, in favor of a canonical ensemble or an ensemble where the pressure 
is the control parameter, however it is difficult to find a convincing argument. It is 
the dynamics of the reaction which determines which fraction of the energy goes into 
collective motion, and which fraction into particle production. This has nothing to 
do with a heat bath nor with constant pressure on the droplet. Consequently, the 
relation between the energy fluctuation, seen in a system with a fixed temperature, 
and the true energy fluctuation is all but evident. 

Therefore, we use another approach to study the influence of energy fluctuations 
on the observables. We assume that the volume of the system remains unchanged in 
order not to have too many variables and that the energy fluctuates. For technical 
reasons, we use discrete distributions, using E% = i AE, with AE = lGeV. For y/s = 
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Figure 8: Average transverse momenta as predicted in the phase space calculation 
as a function of in comparison with the experimental values |18| for different 
particle species. 
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Figure 9: The Gaussian, Poissonian and the NB energy fluctuations. 



200 GeV, we have (E) = 16.15 GeV from the above fitting work, and correspondingly 
we take < i >= 16.15. We study three cases: 

a) The energy distribution is Poissonian 

Prob(z) = (2)teXP , ( - 

and 

b) The energy distribution is Gaussian 



Prob(i) 



^exp(-^) 2 whenEi e [2.5 GeV, oo) 



0.63 v^ct 

otherwise. 



where an energy threshold of 2.5 GeV is taken for the proton-proton system. // = 
a = 14.01GeV to obtain (i) = 16.15 and the factor 0.63 is used to normalize the 
energy distribution, and 

c) The energy distribution is a negative binomial distribution 



> NB (-: ; n,k)- k{k + 1) "' {k + i ~ 1) 



n+k ' 



i\ (n + k) 

The negative binomial distribution is well normalized, and < i >= n. So we take 
n = 16.15. The parameter k = 3 is chosen to get the best fit to multiplicity 
distribution data from UA5. 

All the three types of energy fluctuations are displayed in fig. 9. In Fig. 10, we 
display the influence of these fluctuations on the charged hadron multiplicity distri- 
butions. We compare the results from a fixed energy of 16.15 GeV, the energy with 
a fluctuation of the above-mentioned Poissonian, Gaussian and the NB type. We see 
that already for a fixed energy the fluctuation of the charged particle multiplicity is 
considerable. Different energy fluctuations give different multiplicity distributions. 
The NB energy fluctuation reproduces the UA5 data for non single diffractive events 
in antiproton-proton collisions at = 200GeV. 
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Figure 10: Distribution of the multiplicity of charged hadrons. We display the 
results from a fixed energy of 16.15 GeV, the energy with a distribution of the 
above-mentioned Poissonian, Gaussian and the 2-NB type. Also displayed are the 
UA5 data for non single diffractive events[20j. 
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0.179 


NB E dis. 


0.208 


0.163 


0.241 


0.180 



Table 1: Different K to 7r ratios where (...) means event averaging. 

How does the multiplicity of identified hadrons fluctuate if the droplet energy 
fluctuates? This is studied in fig. 11, where we display the multiplicity distribution 
of the most copiously produced particles for fixed energy E = 16.15 GeV and for 
a Poissonian, Gaussian and NB energy distribution. We see here as well that al- 
ready for a fixed droplet energy the multiplicity fluctuations are important. Though 
different energy fluctuations cause different multiplicity distribution, the energy fluc- 
tuation gives very little effect in the average multiplicities. So our approach, fixing 
microcanonical parameters by fitting the averager multiplicity data, is quite reliable. 

There is also a correlation between the pion and kaon multiplicity in a given 
event, shown in table. 1. The event averaged K/ir ratio is considerably different 
from the ratio of the average kaon and the pion multiplicity. The energy fluctuations 
change more the event averaged K/ir ratio than the ratio of the average kaon and 
the pion multiplicity. 

6 Conclusion 

We have presented a micro-canonical phase space calculation to obtain particle mul- 
tiplicities and average transverse momenta of particles produced in pp collisions as 
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Figure 11: Multiplicity distribution of charged pions and kaons from a fixed energy 
of 16.15 GeV, the energy with a distribution of the above-mentioned Poissonian, 
Gaussian and the NB. The number refers to the mean multiplicity. 
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a function of \fs. Using the multiplicities of p, p, n + , n~ , we fit the two parameters 
of the phase space approach, the volume and the energy. 

Using these two parameters, we calculate the multiplicities of all the other 
hadrons as well as their average transverse momenta. The calculated multiplici- 
ties agree quite well with experiment as far as non strange hadrons are concerned. 

For the yields of strange hadrons (as well as those with hidden strangeness), 
the prediction is off by large factors. In canonical and grand-canonical approaches, 
strangeness suppression factor have been used to solve this problem. 

The energy obtained by this fit is much smaller than the energy available in 
the center of mass system of the pp reaction, because part of the energy goes into 
collective motion in beam direction. Nevertheless, the average transverse momenta 
of the produced particles (not only non-strange but also strange) from this fitting 
agree quite well with experiment. 

We learn that the volume of the pp collision system increases with the collision 
energy. However, it saturates at very high energy(with Antinnuci's parameterization 
as input). The maximum value does not exceed 100 fm 3 . Together with the results 
from article [Hj, we conclude that the grand-canonical treatment cannot describe 
particle production in pp collisions even at high energy. 

We study the effects from energy fluctuations and find that it is quite reliable to 
fix the microcanonical parameters without considering energy fluctuations. 
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